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The problem of determining whether a line of bearing 
measured by a local surface-based sensor coincides with a 
threat whose position has been previously estimated by an 
ocean surveillance sensor is addressed.  Uncertainties in 
the position estimate of the threat, in the bearing error 
and in the position estimate of the sensor are considered 
in measuring the probability that the threat lies on a given 
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I. @ENTRODUCEION 


As long range surface launched weapons systems continue 
to be introduced into the fleet, the operational commander 
is increasingly faced with the problem of being able to 
launch weapons at targets located beyond the horizon. 
Successful employment of such weapons is dependent not only 
on the ability to detect, classify and localize targets at 
considerable distances, but also on the ability to distinguish 
the true target from a potentially larger field of false 
targets. While long range ocean surveillance sensors may be 
of assistance in the identification and localization of tar- 
gets, the information provided may not be refined sufficiently 
to permit effective targeting of long range weapons on that 
basis alone. The on-scene commander must in general rely on 
additional data on target location gathered locally and close 
to the time of weapons launch for accurate targeting. Thus, 
he must still be able to detect and track the desired target 
and be able to distinguish it from other targets within 
range of his sensors. 

The procedures developed in this paper are designed to be 
of assistance in addressing the last of these problems. They 
are applicable when the information available is an error 
ellipse around a threat location estimated by an ocean sur- 
veillance sensor and bearings only data generated by a local 


surface-based sensor. The question of whether a target 


a 





detected by the local sensor is the same as that whose 
estimated position was provided by an external sensor can 
only be addressed if information is available on the loca- 
tions and tracks of all possible targets within range of the 
local sensor. Since such data is generally not available, 
this paper does not attempt to answer that question, but 
rather develops a method by which bearing information from 
different sensors can be compared as to the likelihood of 
each bearing being associated with the threat identified 
previously. It is envisioned that these likelihoods can 

be then used to induce an ordering among bearing data gathered 
by different sensors or, conceivably, conflicting data 
gathered by one sensor. The ordering would be based on the 
likelihood that each bearing will contribute to refining the 
Original estimate of the location of the target of interest. 
This information could be applied in a number of ways: as 

a guide to allocation of more capable sensor resources for 
purposes of obtaining targeting information; as a guide for 
allocation of weapons against more than one threat; as a 
means of pre-processing data before entering it into a target 
motion model, thereby reducing the chance of introducing 
unrelated data. 

To determine the likelihood that a given line of bearing 
and the threat coincide, consideration was given to the 
uncertainties inherent in estimation of target position, 
in the measurement of bearings by a particular sensor and 


in estimation of sensor location. It is assumed that at some 





time to! an Ocean surveillance sensor detects a threat whose 
position is estimated to be within an elliptical region with 
p, x 100% certainty. The estimated position data are received 
and converted by the on-scene commander into a probability 
distribution described by a truncated bivariate normal density 
Function. 

It is further assumed that the standard error 0. charac- 


teristic of the local sensor is known. The sensor bearing 8 


with bearing error c, is then projected from the sensor 


B 
position through the threat density function. 

Since sensor position relative to the target may itself 
be subject to navigation error, the uncertainty is introduced 
as a truncated circular bivariate normal distribution centered 
at location (u Yo? with standard deviation c. 

A TI-59 calculator program is developed which estimates 
the likelihood that the threat identified by an external 
sensor lies along bearing 8, given the threat distribution, 
the bearing error, and sensor position distribution relative 
to the threat. 

The theoretical basis for this calculation is presented 
in Chapter II. The algorithms used in designing the caicu- 
lator program are described in Chapter III. A program listing 
and verbal flow are provided in Appendix A along with instruc- 
tions for the user. Appendix B contains a verbal flow of a 
program designed for use when considerable time has elapsed 


since the initial estimate of the location of a moving 


target. 





II. THEORETICAL BASES 


The general approach to determining the likelihood that 
a measured line of bearing 8 is the true bearing from the 
sensor to the threat identified and localized by an external 
sensor is discussed in this chapter. Calculations required 
when using threat position information both as initially 
generated by the external sensor at time zero, to’ and as 


distorted to account for an intervening time late, t are 


L 


discussed. Uncertainties in bearing measurement and sensor 
position are included. 

Initially, assume that sensor position is known with 
certainty. Let Br be the true bearing of the threat from the 


sensor. Since threat location is uncertain, is a random 


PT 

quantity with probability density function f (Bm). Let 
5TRUE * 

£ (87 Bp) dB be the probability that the errors in bearing 


measurement are such as to give rise to a bearing on the 
threat in the interval d8 about the observed value 8 when 
the true bearing is E Then the likelihood of observing a 


bearing ß is: 


ከ. J £,(878 


| URL) SS (1) 
all B E E 


T  8$vRUE 


The density £f (8m) is determined by the probability 
“TRUE ^ 


density function of threat location at time t, fy y (Xy; t). 


TO 





Given the probability density function of threat location 
at time t and assuming that bearing errors are normally dis- 
tributed with mean zero, consideration is limited to the 
probability that the observed bearing is the true bearing 
of the threat, given the threat is contained in a planar 
region E and the true bearing Sr is in an interval about ß 
with upper and lower bounds Bu and Br The region E is 
selected to be the minimum area planar region which contains 
the threat with a specified high probability Py * The inter- 
val (Bzr By) is selected so that, for all the possible values 


of 8, contained in the interval, 8 is within an interval of 


T 
specified high probability 523 around Ba The fan LBC: Bu 

is symmetric about 8 so that it represents the minimum area 
region which meets the above criterion. The likelihood of 
observing 8 when the threat is in E and Bm is in (8, ; 8g) is 
determined by integrating expression (1) above, over a region 
defined in the manner of the shaded area of Figure l(a). 

If sensor position is not certain, assume that it is 
distributed in accordance with a circular bivariate normal 
distribution Ey y (Uv) - Again, consideration is limited to 
determining the likelihood of observing a bearing g given 
that the threat is contained in a planar region, E, 85 is in 
LB: Bu) and the sensor is contained in a planar region C. 
The region C is selected as the minimum area planar region 
which contains the sensor with a high probability P4. The 


region of E over which f£'(g) is evaluated expands as the 


shaded region of Figure 1(b). 


11. 








(a) Sensor Position Certain 





(b) Sensor Position Uncertain 


FIGURE 1. REGION OF INTEGRATION 
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Since the measurement errors involved in estimating the 
threat position, the sensor position and the bearing angle 
arise from different measurement procedures, it is reason- 
able to assume that errors are independent. The probability 
densities required to perform the above calculations can be 
estimated as follows. 

At time to an ocean surveillance sensor estimates the 
position of the threat to be located within an elliptical 
area characterized by the parameter set E = {X,Y,9,A,B} with 
confidence P] * 100%. The elements of the parameter set E 
are respectively: X the latitude of the estimated threat 
position, Y the longitude of estimated threat position, 

6 the orientation of the major axis from true North, A the 
length of the semi-major axis, and B the length of the semi- 
Minor axis. Since the measurement errors in determining 
threat position are generally assumed to be normally distributed, 
the ellipse characterized by E represents the minimal area 
p, x 100% confidence region about the mean (X,Y). Treating 
this ellipse as the Py probability region of a bivariate 
normal distribution, a density function for the threat posi- 
tion at time tg can be estimated. For convenience, locate 
the origin of a rectangular coordinate system at the center 
of the ellipse, (X,Y), with positive x-axis located along 
the major-axis of theellipse at a bearing 0 from true North. 
Assume a flat earth in the region of interest. Let e 0. 
The mean of the threat position density fy y(Xıy;0) is then 


the point (0,0). The variances in the X and Y directions can 


T 





be derived from the fact that the region with minimal area 
which contains the threat with probability Py is a k-sigma 


ellipse where k is determined from the relationship: 


| =k*/2 
P[threat located in k-sigma ellipse] = l-e 
Thus, 
2 . 2 

Sy A ES) 
and 

NS 
where 

k = -2 in(1l- pi) f 


X and Y are assumed to be independent. 
If the course and speed of the threat are known with 
certainty to be yj and s respectively, the probability density 


of the threat position at time late t, can be shown to be 


L 
again a bivariate normal with mean (st, cos(8-y) ‚st, sin(8-y)) 
and variances Man Ey 0 The Pi probability region of the 


density at time t, would then be an ellipse congruent to that 


L 
characterized by the set E above but centered at the point 


(st, cos (8-y),st sin(8-y)) (Figure 2) [Ref. 1]. 


L 
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ጋ. 
(st, cos (8-4) st, sin(9=1)) 





FIGURE? TIME LATE ELLIPSE, KNOWN COURSE AND SPEED 


IAS 





In some cases, the motion of a submarine on patrol in a 


large area can be characterized, when t, is large, by an 


L 
expansion of the probability area with time at a rate D. 


The result in such a case is that L y Ge Yi tu) is still 
/ 


bivariate normal with mean (0,0) but with variances oy + Dt, 


2 
and oy + Dt, - 


When the motion of the threat cannot be described in the 
above manner, and the course and speed are not known, but 
assumed to be distributed according to the densities f, (4) 
and fo(s), determining the probability density fx yXıyit,) 
is a considerably more complex problem. Let (x, (s,v) vg (sv) 
be the coordinates of the point at which the threat would 


have to be located at time zero in order to reach the point 


(x ) at time t, if the threat speed were s and course y. 


om E 


Thus, x, (s, v) = x, -st,cos(8-y) and yg (s,v) = y,-St,sin(8e-y) 


(rgure 3). Then, 


e — 360 
fy y00Yit = | NS. £y y X5 (8/9 Y (SW) it] 


£y (4) £s (s) dyds [Ref. 1] 


This density is no longer normal. In the special case when 


A = B, i.e., 2:2: = D y has a uniform distribution over 
the interval (0?,360?) and S is known with certainty, as 
shown in Reference 2, the density changes with time as in 


Figure 4. 


l6 








FIGURE 3. TIME LATE POSITION, COURSE AND SPEED UNCERTAIN 
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FIGURE 4. TIME LATE THREAT DENSITY, SPEED KNOWN 
COURSE UNIFORM (0°,360°) 





The distribution of the bearing error measured by the 
local sensor can be estimated if the standard error of the 


Sensor, oc is known. The bearing error is then assumed to 


g 
be normally distributed with mean zero and variance 0. 
Assume that the estimated sensor location (uo, Vo) is 
accurate to within R nautical miles with p4 * 100% confidence. 
The density fy,y UV) of sensor location can be assumed to be 
EE CIrcular bivariate normal with mean (uo, V9) and variance 
T = የ The value of k is determined from the relation- 

ship: 
| | | ! ke 

P[(u,v) contained in k-sigma region] = l-e ; 
where the probability on the left is P4 in this case. 

Since the evaluation of £'(8) considering sensor position 
density and bearing error density does not generalize to a 
closed form, algorithms are developed in the remainder of this 
paper for estimation of the probability (likelihood) &hat the threat 
identified by an external sensor lies on a line of bearing 
measured by a local sensor given that the threat is in region 


E, the sensor is in region C and the true line of bearing 


lies within (Bj r 8) - 


Ike, 





III. ALGORITHMS 


As indicated in the previous chapter, the variety of 
geometrical situations which can arise depending on the 
location of the sensor relative to the estimated threat loca- 
tion precluded development of a generalized analytical pro- 
cedure. Rather, algorithms are developed in this paper for 
numerically evaluating f'(8). The procedures applicable to 
a bivariate normal threat location density have been imple- 
mented on a TI-59 calculator. Appendix A contains a listing 
of that program. The calculations required when the time 
late threat location density is no longer normal exceeded 
the available program capacity of the TI-59 and therefore 
have not been implemented. A detailed verbal flow is provided 
at Appendix B for future implementation on a larger machine. 
The algorithms used in both situations are described in this 


ehapter. 


A.  BIVARIATE NORMAL THREAT LOCATION DENSITY 

This case includes situations (1) where the time elapsed 
.Since generation of the initial error ellipse by the ocean 
surveillance sensor is negligible, (2) where the motion of 
the threat can be assumed to be random in the manner described 
above, and (3) where course and speed of the threat are 
assumed to be known with certainty. With appropriate modi- 
fications to the input data, all three of these situations 


can be addressed using the program contained in Appendix A. 
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In situation (1) the data entered are the parameters of the 
ellipse as generated by the ocean surveillance sensor. In 
(2), the location and orientation from North of the ellipse 
is the same as originally generated, but the size of the 
ellipse expands at some constant rate of area per unit time 
which must be estimated by the user. This rate D, times the 


elapsed time, t yields the factor which must be added to 


Ts 
the seni-major and semi-minor axes of the original ellipse. 
That is, if the original error ellipse is a p, * 100% Conti- 


dence ellipse, the semi-major and semi-minor axes of the 


diffused ellipse will be input as A' and B' respectively: 


E ln (1-p,))Dt, 
CM ስጮስጮ ln (1-p,))Dt, 


In situation (3) the dimensions and orientation of the time 
late ellipse are the same as those of the ee error 
ellipse, but the center of the ellipse is displaced from 
its original position by the known velocity times elapsed 
time. The updated position of the error ellipse is treated 
as the origin of the rectangular coordinate system for this 
situation and all linear measurements are made relative to 
this system. All angular measurements are made from true 
North. 

Estimation of the likelihood that the threat lies on a 


bearing from the local sensor given that the threat is located 


21 





in a p, x 100% confidence ellipse and the true bearing lies 


within the bounds ) with confidence po x 100% proceeds 


(B, 187 
as follows: 


l. Estimate parameters of the bivariate normal density 


u uoces 2 ሮ 2 
X Y : pre Y 


2. Determine Br and Bp such that an interval of length 


25“. centered on either Br, or 8 


f. y GC Y) of threat location: 


U would contain 8, the 


measured bearing, with probability P3: 2B dx 


Pg 86 
= + ረ 
By B K 205 | 
3. Determine sensor location coordinates relative to 
the origin of the threat ellipse. 


4. Subdivide the angular interval (5178 ) into 2n sub- 


U 
intervals. Each subinterval k intersects the ellipse in a 
Strip with average width Wy which corresponds to AB, - 

5. At the midpoint of each subinterval k determine the 
equation of the line through the sensor position at the true 
bearing By from North. 

6. Let the equation of the line of bearing &, be X ==. 
Then the plane perpendicular to the xy-plane which contains 


this line intersects the bivariate normal threat density in 


a curve whose equation is 


2 2 
zu DES) 
m co g 
ao) = ae 3 E (2) 


found by making the substitution X - = in the density 
Eunetion f. y (XY) - It will prove convenient to expand the 
/ 


right side of (2) as follows [Ref. 31: 
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E c 
: ወ ተጉ ው 
muU S e (3) 
“2 Vo ቱ ፲በ o 2 
Y X 
2 2 
l ey Em er Coy 
c ? mo 2 Í 29 e z Ed +m 7 
RS X E d y ur 7% 
ሃሬ 2 “ሇም 
For computational purposes, assume that the width NL 


of the region of the ellipse cut out by the angular subinter- 
val around By will be nearly constant through the ellipse. 
Let the points of intersection of the line of bearing By 

with the threat ellipse be (X 


and (X Then 


k1' * KL? k2'Yx2) - 
approximate the volume of the normal density over this region 
by the absolute value of the product of the area under the 


curve g between Y and Y and W, (Figure 5). Observe that 


kl k2 k 
the term in brackets in equation (3) above is the density 


function of a univariate normal random variable with mean 


and variance 


The term preceding the brackets in (3) is the slope of the 


line of bearing, m, times the density of a univariate normal 


ን ; 2 2 "2 
random variable with mean zero and variance o. +m Ty 


2:3 








e 
2 


FIGURES 
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evaluated at c. Thus, evaluating the expression 


> 


IM. SJ sy) dy 
"kl 


is equivalent to the computation 


l 
|W, m(———— HL = DA O E 
k 2 i: 
NS S : 
Y x SX X 
where 
eo 5 c me = 
_ Y Y X 
MENS 2 IA m ) 
g tmo S 
Y X Ñ 
for i = 1,2, 9 is the Pl[2<z] when 2 is a standard normal 


random variable, and 9 is the density function of a standard 
normal random variable. 

If the substitution Y = mX+c were made in the density 
fy. y UY) an analogous situation would arise with the limits 
of the integration being X] and ፳፲ - 

7. If the value of the density under the bivariate 
normal curve over the region of the ellipse defined by the 
subinterval AB, is weighted by the probability of observing 
a bearing error (By. 7 8) the result is the probability of 
observing the bearing of the threat as 8 when the true threat 
location is in the segment of the ellipse defined by AB - 


This probability is approximated as follows: 
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> 


y-c , 
L Ñ O m 


Kl 


Since the bearing error is assumed normally distributed with 
mean zero and variance Je the value of SELLE can be 


determined by the expression 


where $ is as above the density function of a standard 
normal random variable. 


3 


f'(8) = i £,(8:8,) £ (8 


a ) dß 
all By TRUE 


E ic 


is approximated by 


NIU ee ren a ee gee rane! 


The value of the sum is determined by repeating steps 
A.6 and A.7 above at the midpoints of each of the 2n sub- 
intervals defined in step A.4 and summing the results of 
each of these calculations. Obviously, the finer the sub- 


division of (8 the more accurate will be the estimate 


T^ Bp! , 
of the likelihood, but also the longer the calculation will 


take. 
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Dee SENSOR POSITION UNCERTAINTY 

The result of the above calculation will be the likeli- 
hood f'(8) that the threat lies along the bearing measured 
by the sensor given that the threat is located within the 
threat ellipse and the true bearing is within the interval 


(8 ) and the sensor is at the position used to perform 


To 
the calculation. We next will introduce additional calcu- 
lations that are required to account for the fact that the 
sensor position is not known with certainty. 

l. If the position of the sensor is estimated as being 
within R nautical miles of (u Vo) y its assumed coordinates 
in the xy-system, with pz * 1003 confidence, estimate the 


parameters of the sensor location density f (uv) wren 


mean zero and variance ot; mean = (uo Yo?” variance, 
= EN In(1-p,)) in threat centered coordinates. 


2. The bearing measured by the sensor is 8 regardless 
of the sensor location. Assume that the area of intersection 
of the angular wedge (8, 8g) and the threat ellipse does not 
change significantly as the sensor position is moved along 
the line of bearing 8. Then the effect of the bivariate 
normal distribution of sensor location can be approximated 
by considering only the univariate normal density along a 
line through uer vo) perpendicular to B. Repeat the calcu- 
lations in steps A.3 through A.8 above with the sensor located 
at each of the three points (us Yo?” (u_+.970cos(8-8-90), 
v¿+.970sin(8-8-90)), and (u_7.970cos (8-8-90) ,v_-970sin(8-8-90)). 


If the line through (UL Yo? perpendicular to 8 is subdivided 


ረ 


symmetrically about (ue such that 1/3 of the univariate 
normal density lies above each subinterval, the three points 
chosen above represent the "center of gravity" of each third 
of the density (Figure 6). 

3. If R is chosen to include a significant proportion 
of the sensor density, i.e., on the order of 2c or greater, 
the probability of the sensor being located in each of the 
Ecce regions is approximately iI. Thus, st P4 is on the 
order Of .86, multiply each result in step B.2 by 1733 

4. Summing the results of steps B.2 and B.3 yields an 
estimate of the likelihood that the threat lies at bearing 
8 given the threat is in the Pp, x 100% confidence ellipse, 
is in (8 


8 ) and the sensor is in a p4 * 100% confidence 


T pus 
region. That is, the likelihood is estimated by 


where j is the index of sensor position in figures 5 and 6. 
Instructions for application of Enhe TI- 59 progranito 
calculate the above are included in Appendix A. 
C. THREAT DISTRIBUTION NOT BIVARIATE NORMAL AFTER TIME 
LATE ELAPSED 
The basic approach to evaluating £'(8) when the time late 
distribution of the threat is no longer bivariate normal 


is the same as that just discussed. The principal difference 


28 





din 


2 
/ 
x m 
R nn 
(u Vo -(u, v) 
eT 
(u,7.97ccos (678-90) ; (u_+.970c0s (8-3-90) ፡ 
ኛዕ”» 272848) (5-2>30)) OT sin(3-8-90)) 


= (u4, V.) 


= (9372) 


FIGURE 6. ESTIMATE OF SENSOR POSITION DENSITY 
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arises from the fact that the time late density is signi- 
ficantly more complex in this case. 

The threat distribution becomes distorted from the normal 
after some time late when the course and speed of the threat 
are constant but not known with certainty. Application of the 
method described herein requires that the user assume a dis- 
crete distribution of the speed of the threat with upper and 
lower bounds Sy and Sy, respectively. In addition, the threat 


course is assumed to be uniformly between 0? and 360?. 


The density of threat location after some time late t 


L 
when the threat speed is S; then becomes [Ref. 1]: 
360 x-s,t, cos (8-4) ) ? 
£ [መሇያ s.) = ms f exp [-3( INE 
pay ttn" Dic. 2 2 
AY O, 
X 
(y-s.t sin(e-4y))* 
mnc A o yy EY 
eZ 360 
Y 


where 9 is the bearing of the major axis from North. Note 
that the new threat density is still centered at the same 
position as the time zero ellipse but its shape changes as 
in Figure 4 of Chapter II. If the threat speed is Sir and 
the course is uniformly distributed ever (0-, 260°), the 
outer limit of the new planar region containing the threat 
after time tr, has elapsed, given that it was originally 
located in the Pi * 100% confidence ellipse with semi-axes A 
and B, can be represented by an ellipse with semi-major axis 
(Figure 7). Thus the 


Ats,t and semi-minor axis B+s.t 


L L 
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ወ mE 





mE 





region over which the density will be evaluated is still 

the intersection of an ellipse with an angular wedge, 

recalling however that the density is no longer normal. 

Further, the new elliptical region does not represent a 

Pp, x 100% confidence region of the time late density. 
Calculation of the likelihood that the threat lies along 

bearing 8 given that the threat lies in the 51 * 4005 ellipse 

at time zero, that the true bearing lies in (Bir By)? that 

the sensor lies in the P3 * 100% circle and that the speed is 

S; proceeds as follows: 


l. Estimate the parameters of the original normal dis- 


tribution: mean = (0,0), 2." = ። eee In(l-p,)), 
2 . 
E B^/(-2 ln(l p,)). 
2. Determine the upper and lower bounds on the true 


bearing wedge, 8 +k,9, and 8 7 Ka0g- 


3. Determine the time late planar region as the ellipse 


with semi-major axis equal to A+s,t, and semi-minor axis 


equal to B+sj¡t,- 
4. Determine the position of the sensor relative to the 
ellipse center. 


5. Subdivide the bearing fan (8=k,0¿1/8+k,0,) into 


2n subintervals. 
6. At the midpoint of each subinterval, determine the 
equation of the line through the ship position at that 


bearing By 
7. Let the equation of the line of bearing By be 
x = ፍር. Then, the plane perpendicular to the xy-plane which 
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contains this line intersects the time late threat density 


in a curve whose equation is 


Er 
2 2 2 2 2 
g ENC g Taille © 
g, (y) m E Y X Mi እሪ X 
57 መ + 829." vau m Sy Oy 
1. %y +m ሠ 2 c s 2 
1 m “ሄ oy Oy SS Ox 
360 (s,t,) cos” (9-4)-2s,t, ycos (9-4) 
E A 
X 
ጋ ጋ 
8 (sitz) sin Er ae ame E N 
2 360 


Cy 


found by making the substitution x = <= in the density 
Ennetion fy y Yit) and expanding. 

For computational purposes, assume that the width Wi 
of the region of the time late threat ellipse cut out by the 
angular subinterval around By will ee nearly constant through 
the ellipse. Let the points of intersection of the line of 
bearing Br with the time late threat ellipse be (Xer Yk)” 
(Xar Ypa). Then approximate the time late density at speed 


s, over this region by the absolute value of the product of 


W, and the area under the curve gr between Yel and Yo” The 


k 
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area under curve g, can be approximated as follows.  Sub- 
divide the interval (Y, Ta) into nj segments of Length h. 
Evaluate 9, at the midpoint of each segment, Y Note that 
the first term in parenthesis in (4) is a constant equal to 
m times the value of the density of a normal random variable 
with mean zero and variance oy mo," evaluated at c. The 


second term in parenthesis is the density function of a normal 


random variable with mean 


and variance 


evaluated at y. The variable y also appears in the integral 
term in equation (4). Numerically evaluate this term of 
(4) with y = Yr. Let 91 (Y y) be the result of this computation. 


Then, evaluating the area under the curve Sr between Yel 


and Ye is equivalent to the calculation: 
p - d o, ^ * mo, 
T ĖS 
DA Ce d 
F X Y X 
where ¿(-) is the density function of a standard normal random 


variable, and 
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An analagous situation arises if the substitution Y = mX+c 
is made for y in the time late density fy yX Yit). 

8. The value of the density over the subinterval con- 
taining By is then weighted by the instantaneous probability 
that the bearing error is By 7 8: 


7x2 


Wy p TEE EE 
kl 


9. Steps C.7 and C.8 are then repeated for each sub- 
interval of (87 8p) and the results of each calculation are 
summed. 

10. The uncertainty in sensor location is accounted for 
by repeating steps C.4 through C.9 with the assumption the 
sensor is located at each of the three points in Figure 5, 
multiplying by the probability the sensor lies in that 
interval and summing each result. 

11. The result of calculations in steps C.l through 
C.l0 is the likelihood that the threat lies on bearing ß 
given the threat was originally located in the p, x 100% 
confidence ellipse, E € (8p, Bg), the sensor is located in 
the pP, x 100% circle and the threat speed is Si- The con- 


dition that the speed is Si is removed by repeating the 
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calculations C.l through C.10 for each of the speeds Sir 
Do-c1,...,M multiplying thegdnsesubtabthe probability that 
the speed equals Si and summing all M results. The final 
result is the likelihood that the threat lies on 8 given 

it was originally located in the p, x 100% ellipse, 

Bp € (8r, Bg), and the sensor is in the Pp, x 100% circle esti- 


mated as 
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IV. TCONCCUSTONS 


Possible applications of and extensions to the algorithms 


developed in Chapter III are discussed in this chapter. 

As indicated in Chapter I, the objective of this paper 
has been to develop a means of assessing the likelihood a 
threat whose position has previously been estimated lies at 
a given measured bearing from a local sensor. The procedures 
were developed with a view towards permitting the user to 
make comparisons among lines of bearing measured by different 
sensors or among conflicting bearing information generated 
by one sensor. The approach chosen has been to estimate the 
likelihood that a threat lies on bearing 8 given that the 
threat is located in an ellipse of specified confidence 
pi * 100%, that g is measured with Po * 100% accuracy, and that 
the sensor position is measured with Pz * 100% accuracy. The 
algorithm to calculate the likelihood in the cases where the 
probability distribution of the target can be assumed to be 
bivariate normal at the time of the bearing measurement has 
been implemented on the TI-59 calculator. The cases in which 
this program applies are the following: (1) when the time 
elapsed since generation of the threat error ellipse is small 
enough to justify using the original estimate of the ellipse; 
(2) when threat course and speed are known, in which case 
the ellipse center is translated from the original position 


according to the course, speed and elapsed time information; 
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(3) when the threat can be assumed to be moving about in a 
random manner over a significant region in such a way that 
the ellipse center remains unchanged, but the x and y vari- 
ances have increased. When none of these cases hold, but 
the course is assumed uniformly distributed over (0°,360°) 
and speed has a discrete distribution over a finite interval, 
the threat density at the time late tr is not a bivariate 
normal density. The algorithm applicable in this case has 
not been implemented, but is described in some detail in 
Chapter III and Appendix B. 

Once the appropriate computation has been completed for 
each of the bearings considered, the results can be used to 
weight the value of several bearings in refining the threat 
location estimate provided by the external sensor. Note that, 
although unlikely, £'(ß) may correctly be greater than one. 
Comparisons using these likelihoods should be made only when 
the upper and lower bounds on the true bearing fan for each 
bearing are chosen at the same probability level P2 and the 
uncertainty areas for all sensors include the same probability 
level Pz: Further, these likelihood levels should be selected 
so as not to exclude a significant portion of the appropriate 
density. ፲፻ pj Or p, are not the same in all cases to be 
compared, f'(8) must be divided by the applicable value of 
Po Or P3 for each measured bearing 8 to be considered. 

Having established the relative value of available 
bearing information, the user can allocate weapons or further 


search effort accordingly. However, the probabilities calculated 
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are strictly ordinal data and do not define a redistribution 
of target location probability based on additional informa- 
tion. Further, the threat ellipse does not contain the 
target with certainty. The power to predict the probability 
of success of a search or weapons allocation plan based on 
the priorities established by these procedures is limited 

by these constraints. In this area in particular further 
research would be useful. 

In situations such as that for which the procedures in 
this paper have been developed, where a track has not been 
developed on the target, introduction of unrelated bearing 
data to a target motion model could impact significantly on 
the reliability of future position predictions. If there is 
high confidence in the reliability of the estimate of the 
threat ellipse provided by the ocean surveillance sensor, 
the prioritization established herein could be used to pro- 
cess bearing data prior to input to a target motion model. 
Using a pre-established threshold, only those bearings which 
coincide with the threat with an acceptable level of like- 
lihood could be used to refine or update a track on the 
threat. 

Desirable enhancements to the algorithms include providing 
for the instances in which the interval of uncertainty of 
the target course is known to be less than (0?,360?). Further, 
if the circular region of radius R contains the sensor with 


significantly less than 86$ confidence or the assumption of a 


39 





bivariate normal distribution of sensor location is unsatis- 
factory, it is left to the user to modify the calculations 
accordingly. 

The utility of these algorithms would also be improved 
by implementation on a larger and faster System than the 


TI-59 calculator. 
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APPENDIX A. 


part I 
Step Number 


woo = 029 
020° 029 
0802-052 
022058 
059063 

064 - 073 
074 - 089 


TI-59 PROGRAM VERBAL FLOW AND USER'S INSTRUCTIONS 


Verbal Flow 


Enter the confidence level P1 EOL 
the threat ellipse. Calculate the 
value of k for the given pi: 


k = v=-2 1n(l-pi) f 


Enter the length of the semi-major 
axis, A. Calculate oy = A/k. 
Enter the length of the semi-major 
axis, B. Calculate o,, = B/k. 

2 
Calculate B/A, o Tyr Tx and Oy 
Enter orientation of semi-major 


axis, @. 


Enter bearing from sensor to center 
of threat ellipse, a. Calculate 


Dey 


Enter distance r from sensor 
position to center or threat 
ellipse. Determine rectangular 
coordinates of sensor position 
(U,V) From polar coordinates 


(-r,9-a). 


Store the constants 360, 2m . 


Initialize register 35 to 0. 
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00° = 096 


እከ 7, = 104 
ከን ን - 143 
144 - 146 
I. X51 
Dl 


Enter number of standard devia- 
tions desired for bearing fan, ka. 
Enter standard deviation of bearing 
error, "5 Calculate Kg 

Enter the angular stepsize desired, 
48, for incrementally stepping 
thronus (8c By) - Calculate žag 
Calculate the largest number n of 
increments of size A8 contained in 
Ko. - 558 degrees.  Initialize 
counter 00 to n+l. Save n+l in 


register 20. Determine 


E rl T 
SAB = K Og 558 lis 


the residual increment. Calculate 


T 
25948 


Initialize counter 0I to 2. 
Calculations will be made at the 
midpoint of each interval from 


o NO 


ese of ee from 

seo TAB o) c Koog , and finally 
at 8. Counter 0l indicates whether 
calculations are complete on both 


sides of 8. 


Enter bearing measured by sensor, 
Be 


Enter standard deviation of sensor 


ገርን 1.11 31፡72 
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155 


mo - 157 
Mo = 238 
ee = 259 
mow = 340 


Enter index of the sensor position 


to be used for this run. 


Coordinates of the sensor position 
are selected in accordance with 
run number entered above in 


Subroutine sin. 


Determine whether a bearing 8' 
parallel to either axis is included 
in the fan (8; (Bg) - Tr (8, , 8g) 
includes a bearing parallel to the 
y=aáaxis TUSE program Irori Parco 
The appropriate program number is 
displayed in calculator display 


register. If (51፡8 ) does not 


U 
contain a bearing parallel to 


either axis, use program l. 
Subroutine P/R. 


Subroutine sin. 
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Part II, Programs 1 and 2 


Step Number 


MOD — 003 
004 - 007 
Dur = 012 
ሸም» 020 
ከር ን 028 
IDE 050 
031 - 045 
Udo = 051 


Verbal Flow 
Initiate eg’. 


If the last angular increment on 
this side of 8 has been considered, 
go to step 513. Otherwise 


continue. 


Decrement counter for angular 
increment. rf counter - 0 ር) 
021. Otherwise continue. 


Remove flag to indicate this is 
not the last angular increment. 
Recall 48, the input angular step- 


size. Go to 029. 


Set flag to indicate this is the 
last angular increment on this 
side. Add one-half the input 
angular stepsize Ag and one-half 
the residual stepsize 828: 


1 H 
258 E 2968» 


Increment 8' by the appropriate 


stepsize. 


Convert 8' to an angle between 0° 
and 360°. 


Calculate 9-8'. Print 8-8’. 
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052 — 064 


065 - 074 
025° = 086 
ES7 - 088 
EH — 111 
ከ. ን = 132 
MOS 138 


Le =e" |p ae ቢውን ቢ፡፡ ።ዣ  0ዐ።”ን›) 236. 


፲፻ |6-ዩ! | = 0° or = 180°, go to 075. 
Otherwise go to 133. 


If the absolute value of the y 
coordinate of sensor position is 
greater than the length of the 
semi-minor axis, 8, go to 004. In 
this case g' does not intersect the 


error ellipse. Otherwise continue. 


Set flag 2 to indicate that the 
bearing 8' is parallel to the 


x-axis. 


Calculate the coordinates of the 
points of intersection of 38' with 
the threat ellipse. The y-coordi- 
nates are equal to V, the 
y-coordinate of sensor position. 
x-coordinates are determined in 
Subroutine v. The points of 
intersection are symmetric about 
the y-axis. Thus, X] = -X 2" 
Store the smaller x value in 
register 27, the larger in register 
Zon 


Save locations of Xe X27 YLL 


and. TT 


S Y 


Remove flacs 2 and 3 to indicate 


that 8' is not parallel to either 


the x7 or y axes: 
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59 - 


146 - 


55 = 


1 -= 


A ./፦ 


145 


154 


176 


2323 


247 


Calculate the slope of ß': 
m = tan (9-8'). 


Calculate the y-intercept of 8': 
c = Voces. 

I c? E A ^m^ us Se go to 004. In 

this case, 8' does not intersect 

the error ellipse. Otherwise 

continue. 

Calculate X1 and X2! the x- 

coordinates of the points of inter- 


section of g' with the ellipse: 
B c/u 2 2 2 


-mc + = Am +B -c 
X = A 
kl Am? ገ p? 
A? 
-mc ~ — Aam‘ 4 B? - c* 
Pd Iz 
k2 Am? p B? 
A* 


Calculate the y-coordinates of the 
points of intersection of $'with 


the ellipse, ን. Yo 


Led = mX 1 t C 
ኝፔጋ = mi, + C 
Go to 299. 


If |U|, the absolute value of the 
x-coordinate of sensor position, 
is greater than the length of the 


semi-major axis, A, go to 004. 8' 
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248 - 


E50 - 


n 


299 - 


249 


22 


298 


228 


does not intersect the error 
ellipse in this case. Otherwise 
continue: 


Set flag 2 to indicate that Bi is 
parallel to the y-axis. 


Calculate the coordinates of the 
points of intersection of B with 
the error ellipse. The x- 
coordinates are equal to U, the 
x-coordinate of sensor position. 
The y-coordinates are determined 
in Subroutine LONE The points of 
intersection are symmetric about 


መፍ 
==> 


the x-axis. Thus, FE Yun» 
Store the smaller y value in 

register 29, the larger in 30. 
Save the locations of Y Y 


X 


ከን ን kii 


ki? yt 


Determine the width of the strip 


around 8', W: 


If flag 0 set, 


a 1 
W = (d,+d,) 588 (25 AB), 
Otherwise 
a 1 
W= (4,+d,) tan 7AB) . 


di and d. are the distances from 


(U,V), the sensor position, to the 
intersection points (X17 Y) and 


oz 
are calculated in Subroutine log. 


), respectively. Distances 
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e29 = 334 


When using program 1: 


B5 - 380 
E 391 
992 - 401 


Multiply W by A8 if flag 0 not 
set. Otherwise, multiply by $48. 
Save result in register 37. 


IE Br is parallel to either the 
x- or the y-axis, go to 468. 


Otherwise continue. 


Express the y-coordinates of the 
intersection points as standard 


normal random variables, 21 and Z5: 


Sort the values 24 in descending 
order. Let 22 be the larger 
value. 21. is the smaller. 


Calculate $(25') ፦ 9(2,%), the 
probability that a standard normal 
random variate lies between 21 | 


፣ 
and 22 $ 


48 





402 - 424 Multiply the results of steps 
O O 


B )x|m| x W, 


O +M op 
where ó is the standard normal 


density function. 


The result of the calculations in steps 335-424 is 
X 


k2 y-c 
W | Ey v La NY dy 
Yk1 
When using program 2: 
BDD = 372 Express the x-coordinates of the 


intersection points as standard 


normal random variables Z1 dnce ze 


2 
mco | 
ES: 2 
k2 
0 ie m 7 
Y ex 
LEE ; 
» gg 
PE 
Y 2 2 
oO, + Mo, 
i= 1,2 
B73 - 383 Sort the values of 24 in descending 
order. Let 22 = larger value. 
Z,' = smaller. 


il 
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Bear 393 Calculate $(2,') - 9(2,'), the 
Probability that a standard normal 


random variate lies between Z.' 


I 
and 22 - 
394 - 413 Multiply the results of steps 
384-393 by: 
: 5 - ) x W, 
2 2 2 Zu? 
oy + m Jy Sy ግ] Ox 


where $ is the standard normal 


density function. 


The result of the calculations in steps 335-413 is 


> 


W = Ey y (Esmx+e) dx 
kl 


Regardless of which program is in use: 


425 - 465 Multiply results cf previous 
calculation by 


Mi g"-g 
— 6$( ) 
E "5 
" Z 6! 7 65 /< Teor 
MS Sees =260° , 6° > 180% 


and ¢ is the standard normal density. 


The result of this calculation is 


ኃ y 
£,(8:8') |W f £y y CY) dy | 
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466 - 467 
259 — 507 
Eu 512 

SA OS 
316, =2520 


Go to 508 


Calculate the area under the curve 


formed by the intersection of 8' 


with fy y Ue. When 8' parallels 
r 

the x-axis this calculation 

becomes: 

1 


V 
RICO ALIC - ፥(2))] . 


where Z, = Xp o/ oy and 51 = X 1/9 7 
V is the y-coordinate of the sensor 
position. 

When 8' is parallel to the y-axis, 


the calculation is: 


1 U 
where Z and Z 


ን ከር ey M 
U is the x-coordinate of the sensor 
position. 

Go to 422 to complete calculation 
EE OSEE DW [£y 4 06 Y) | 


B 
(GO to 4100'in Program 2). 


Accumulate the probability at each 


angular interval. Display result. 


If flag 49 መር gO to 558. 


Otherwise continue. 


If flag O is not set, that is if 
the angular increment just con- 
sidered was not the last on this 
side of B, go to 004 and continue 
calculation on same side of 8. 


Otherwise continue. 
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521 - 524 Decrement the counter Ol. If 
counter is now 0 go to 543. In 
this case, the probabilities have 
been calculated on both sides of 8. 
The calculation at 8 remains to be 


done. Otherwise continue. 


tig, 


325 = 542 Remove flag 0. Multiply 48, 2 


and =6 AB by -l. Reinitialize 
counter 00 to n+l. Go to 000 to 
begin calculation on second side 
OL Mac 


543 - 567 Set flag 0l indicating that calcu- 
lations on both sides of 8 have 
been completed. The next iteration 
will do the calculation at 8' = 8. 
Remove flag 00. Initialize g' to 0. 
Recall 8. Go to 029. 


558 > 566 Remove flag 01. Display the 
accumulated likelihood. STOP 


The result of this calculation is approximately f, (8) with 


sensor at (U,V) where j = run number entered in Part I. 


567 =- 587 Subroutine y calculates the points 
of intersection when g' is parallel 


to either the x or y axis: 


es os 
ርን END 
X] = -A 71 - (5) 
_ EU 
X2 = ሕግህ (፪) 
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Parallel to y-axis: 


wae CC 
Ee A - D? 
Yu 7 Byl- (^ 
589 - 610 Subroutine log calculates the 


distance from sensor position to 
Pone Or incersection of g' with 
the ellipse: 


ae oe 2 M 2 
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Part III 
Step Number Verbal Flow 


000 - 012 Enter result of run Ll. 
Multaply by 1/3. 


013 - 021 Enter result of run 2. 
Moltiply by 1/3. 


022 - 029 Enter result of run 3. 
MISEEIDIVSSv 1/3. 


0S0 = 032 Display likelihood. 


033 - 041 TE Po is the same for all bearings 
to be compared, enter l. Go to 
042. Otherwise, enter Po for this 
bearing. Divide likelihood by Po: 
Display result. 


042 - 050 I P> is the same for all bearings 
to be compared, enter l and STOP. 
Otherwise, enter P3 for this bearing. 
Divide likelihood by P3- Display 
result, STOP. 
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USER'S INSTRUCTIONS 


The program to determine the likelihood that the threat 
lies along bearing 8 given the threat is in the confidence 
ellipse, Bn E (Bir By? and the sensor is within the P, x 100% 
confidence circle is in three parts. All parts require the 
use of a printer and the use of the Applied Statistics 
Library Module. Prior to running, the calculator must be 
repartitioned: 

l. Enter 4 


2. Press 2nd OP 17 


Bare I 

l. Read sides l and 2 of Part I 

2. Read side 4 of Part II either program 1l or program 
2. Since program 2 is used more often, unless it is known 


that the bearing fan ( ) contains a bearing parallel to 


Bn Bg 
the major axis, recommend using program 2 of Part II. 

3. Press RST 

4. Enter pir the confidence level of the threat ellipse. 

Press A 

5. Enter A, length of the semi-major axis. If data 
provided is length of the entire major axis, divide by 2 
before entering. 
FresszB, 

6. Enter B, length of the semi-minor axis. If data 


provided is length of entire minor axis, divide by 2 before 


entering. Press C. 
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7. Enter 6, the bearing of major axis from North. 
IS ፡ ሜ ቢመቁመዌመሜመሚመጠሜጠመጮልስዶ ሞስ tomentered im degrees. Press D. 

8. Enter a, the bearing of the threat ellipse center 
meen) ener estimated Sensor position. 0° <a < 360°, in 
degrees. Press E. 

9. Enter r, the range from the estimated sensor posi- 
tion to the center of the threat ellipse. Press 2nd A'. 

10. Enter K the number of standard deviations desired 
in one direction from g. The program will construct a fan 
of equal size on the other side of 8. Press 2nd B'. 


ll. Enter c the standard deviation of bearing error. 


a? 
Press 2nd C'. 


12. Enter ag, the desired angular stepsize. Press 2nd 


p. 

13. Enter 8, the measured bearing. Press 2nd E'. 

14. Enter o, the standard deviation of sensor position. 
Press R/S. 


i Enter the number of this run, 1, 2 or 3. When run 
number = 1, sensor is located at its estimated position 
(uo; Vo). when run number = 2, the location will be 
(ut .97ecos(6-8-90),v t .97csin(e-8-90)). When run number = 3, 
1 ።ሪ:፡ን፡ህ፡:፡ቡ!እጅጃጆጂያፕዯያታ፡ል ሰሮርሮ ምን ጋን ርረጋ::ርል ስከ) ve-:975sin(e-g-90)). 
Press R/S. If program number displayed matches that of the 
Part II side 4 read in, continue to 16. Otherwise, press 
2nd CMS, RST. Read in side 4 of the Part II program which 


corresponds to the number displayed. Repeat 3 through 15. 
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16. Press 4 2nd WRITE.  Rerecord side 4 of the Part II 
program read in. This enables data entered in Part I to 


be transferred to Part II. 


Bart Il 

l. Read sides 1, 2, 3 and 4 of Part II Program l or 2 
as selected by the Part I program. 

2. Press RST 

3. Press R/S. Values of 8, and ደክ) 
will print alternately. Final result also prints out at end 
Ot calculations 

4. Record final result: Likelihood threat lies along 
bearing 8 given threat is in p, x 100$ ellipse, true bearing 
js in pP, x 100% fan and sensor is at location used in this run. 
Parts I and II must be completed 3 times (Run numbers 1l, 2 
and 3) before proceeding to Part III if uncertainty in sensor 


position is being considered. 


Bare LI 
l. Read side 1 of Part III. 
2. Enter result of run l. Press A. 
3. Enter result of run 2. Press B. 
d XEUEpntereresubtooborun 3. Eness C: 


ጋ. ስ ከ ሰሜም P5; confidence level of (8,78 ), if necessary. 


U 
Otherwise enter 1l. Press D. 
6. Enter Pa confidence level of sensor position, if 


necessary. Otherwise, enter 1. Press E. 
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APPENDIX B. TIME LATE VERBAL FLOW 
VERBAL FLOW 


1 Pi: Determine the value of k: k = /72 1n (1- pi) 


2. Enter A. Determine m = (እ/ጄ) * 


3. Enter B. Determine Oy 


4. Enter 0, orientation of major axis 


(B/k) * 


5. Enter a, bearing from sensor to ellipse center 
6. Enter range r from sensor to ellipse center 


7. Calculate (uo, Vo) the coordinates of mean sensor 


position. u,= -r cos(0-a), v4- -r sin(ée-a) 


8. Enter Oa: the standard deviation of bearing error 
9. Enter Kae the number of standard deviations to be 
included in the bearing fan on each side of a. 


10. Enter desired angular stepsize, Ag 


ll. Calculate the number of iterations of size Ag required 
_i 
B 2 


the greatest integer less than or equal to the value within. 


on each side of 8: I= [(፳.። A8)/A8] where [-] means 


302" Zu8)/28 is not an integer. Determine 


the size of the fractional increment: 


12. In general (k 


Sg (ikog - 348)/A8 - 1) A8. 
13. Enter g, the measured bearing. Let PROB2=0, PROB1=0 
14. Enter the time late tr 
l5. For each of the discrete threat speeds to be con- 

sidered, repeat steps 16 to 43. Then go to 44. 


l6. Enter the target speed s. Let PROB - O0 
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lj se expand tne ጋስ ን limit of the threat ellipse: 
Let A! = Ar st, 


: = 
B B+st, 


18. Let 8° = 8 
ioe For ፎር ወጋ the T ancrements of size Ag repeat 


stepsez0 to 38. 


moe Let 8° = 8° +46 
ae. calculate (O at). If |e=8" | = 0% or 180°, go to 27. 
we B'| = 90°Zer 270°, go to 30 


22. Calculate m = tan (9-8'), the slope of the line of 
bearing. 

23. Calculate c = mU+V, where (U,V) is the sensor loca- 
tion for this iteration. c is the y-intercept of the line 
of bearing. 

24. Calculate nana p ac If this quantity is less 
than zero, the line of bearing g' does not intersect the 
threat ellipse. Go to step 20 and process next increment of 
size A8, if any remain. If all I intervals of size A3 on 
this side of 8 have been considered go to step 40. 


25. Calculate the points of intersection of the line 


of bearing with the expanded ellipse: 


(-mc 4 (B'/A')VA' ^m? * B'^ - c^) (a!) 


= z A' ^m? t Bi? 
(-mc - (5"/እ')#ሃል! +B? = a) (A' ^) 
A ማያ oo a - 
ጅጋ እ ቪክ” + ቺ!" 


59 





26. GO 32 

27. If |V| > B', the bearing 8' parallel to the x-axis 
does not intersect the threat ellipse. Go to step 20 and 
process next increment of size AB if any remain. If all I 
intervals if size A8 on this side of 8 have been considered, 
go to 40. 


28. Calculate the points of intersection 


eS 
ዘ 


i -B'/1 - (v*/B"*) 


S o 


ጋሩ 
| 


ሪ ን ን ሠ 069). 

30. If |U| > A', the bearing 8' parallel to the y-axis 
does not intersect the threat ellipse. Go to step 20 and 
process the next increment of size 48 if any remain. If 
all I intervals of size Ag on this side of g have been con- 
sidered, go to 40. 


31. Calculate the points of intersection: 
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32. Calculate the median width of the strip of the 
ellipse defined by the angular subinterval under considera- 
tion. If the subinterval is of size Af: 

d +d, 


_ i 
W = 5 (2 tan 255) 


If the subinterval is of size 648: 


d. +d 
W = 2 tan 5548) 


In these expressions d; is the distance from the sensor to 

che ith point Of intersection, i = 1, 2: di = V (0-X,)* * (-,)* 
BLE ።ጭ።ቱ ፡፻ኘ፦፡ 0. በመመ 26. 12 |9-6"| = 90* 

oe 2/0, gO tO 35. 
34. Evaluate the target density from ኛ1 | tO Y,’ along 


the line y = mx + c: 


Let Yi = min (Y, ,Y,) 
i = 
Y, = max (Y, ,Y,) 
ቆ(2) = density function of a standard normal 


random variable evaluated at z 
= 2 2 2 2 2:72 
፳ = (m/ yoy +m Oy ) t(c/yoy +m Oy ) 


Subdivide the interval (YAT TAT) into n segments of length 


h. At the midpoint of each segement, Y; compute: 
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ጋ Soy TA 
ከጀ ( Oy +m Ore / Wo. o.) "EE cona Oy em Iy / ma. o.) ) 
oy +m 9% 


360 1 (st, ) 2 os^ (6-y) - 2st, y. cos (8-1) 
exp [- 5 AA >> 
J L z | ሪጋ 
። m dy 
2 AUAM. 
(st.) sin (6-y) -2st,y. sin(e-y) 
f d V. y Sy 
360 
“y 


when the integral must be numerically evaluated. 
Sum the results over all n segments. 
Go to 37. 
35. Evaluate the target density from Y," to Y.' along 


| 2 
the line x = U: 


Let Yi' = min(Y;,Y,) 
Yo! = max (Y, ,Y,) 
¢(z) = density function of standard normal random 
variable 
Wm i. 
AY 


Subdivide the interval (Yi ', Y!) into n segments of length 


h. At the midpoint of each segment, EE compute: 


360 2 
T (st.) cos(e-y) - 2st. Ucos (e-) 
Rees f el A AAA 
LAO 0 g 
X 
Zu 2 : 
P (st,) sin nn). ay 
360 
EY 
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where the integral must be numerically evaluated. 
Sum the results over all n segments. Go to 37. 
36. Evaluate the target density from X E T 


2 
the line y = V: 


—— 
Let Xi - min(X,,X.) 
X, = max(X,,x,) 
¢(z) = density of standard normal random 
variable 
፪ = = (2 
Y 7 


Subdivide the interval (X,',X,') into n segments of length 


li 2 
h. At the midpoint of each segment, 54 compute: 


2 
X. 360 (st.) cos(0-y) - 2st. x. (cos (8-) 
me ELE | EM Xx 
ox x ወጊ. 
. (st; ) “18 (6>ሀ) ሾር. dy 
2 360 
Y 


where the integral must be numerically evaluated. Sum 
results over all n segments. 

37. Multiply the value of the target density just 
computed by W. 


38. Multiply this result by 


il 


O O 





ፒጩ 
ሯጩ 


where $ is the standard normal density function. 
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39. Let PROB = PROB+ (the results of the calculations 
in steps 20 through 38 for each if the I subintervals of 
size AB). 

40. If the fractional interval of size saAg on this 
side has been considered go to 41. Otherwise, let 
B' = 2668 + Lag. Repeat calculations 21 through 38 once. 
Let PROB = PROB+ (the result calculated at this step). 


41. If the computations on both sides of g have been 


computed, go to 42. Otherwise repeat the computations 


from 18 to 40 on the other side of 8 by letting Ag = -A8, 
SAB = ሙርሰላ. 
42. Let g' = 8. Repeat steps 21 through 38 once. 


Let PROB = PROB+ (the results of this calculation). 
43. Let j = the number of discrete threat speeds to 
be considered. Let PROBl = PROBI + P[S=s]PROB. Go to 15. 
44. Repeat steps 15 through 43 once with the sensor 


located at each of three points: 


(1) (U,V) = (u Vo), the mean of the sensor 
density; 
(2) (U, V) = (u, + .970cos (9-8-90), 


vg * -97csin(e-8-90)); 


(3) (U,V) 


ሮ 


.97 cos(0-8-90), 

v_-.97 sin(8-8-90)); 
where c is the standard deviation of the sensor density. 
The value of PROBl calculated at each iteration will be 


weighted by the approximate probability that the sensor is 
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located in the region of the sensor error circle repre- 
sented by the applicable value (U,V). ፲፻ P3 is of the 
order .86 or greater, multiply by 1/3. 

45. The value of PROB2 calculated after completion 
of step 44 is the relative probability that the threat 


lies along bearing 8 at time t given the threat was in 


Tod 
p, * 100$ ellipse at time tg; B, is in LB: Bu) and the 


sensor is in a p4 * 100$ confidence circle. 
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CALCULATOR PROGRAM 
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